Monte Carlo determination of the phase diagram of the double-exchange model 
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We study the phase diagram of the double exchange model, with antiferromagnetic interactions, in 
a cubic lattice both at zero and at finite temperature. There is a rich variety of magnetic phases, 
combined with regions where phase separation takes place. We identify phases, intrinsic to the cubic 
lattice, which are stable for realistic values of the interactions and dopings. Some of these phases 
break chiral symmetry, leading to unusual features. 
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I. INTRODUCTION. 

The main magnetic properties of doped manganites 
and other materials are described in terms of double 
exchange interactions []4^|. An extensive body of work 
has been devoted to the understanding of this model, 
and the related ferromagnetic Kondo lattice, since the 
pioneering work of de Gennes (for a recent review, see 

0). 

While the model is extremely simple to formulate, as 
it can be written in terms of two dimensionless parame- 
ters only (see next section), the elucidation of its phase 
diagram is yet to be completed. Initial studies || sug- 
gested that the phase diagram could be understood in 
terms of a ferromagnetic and an antiferromagnetic phase 
with an intermediate phase with a canted arrangement 
of moments. More recent work showed that in the region 
of the phase diagram where canting was expected, phase 
separation is more likely |§Hl^l- By now, there are ex- 
tensive studies of the phase diagram of the model in one 
dimension P,|Tl ,|l3| 
dimensions |ij,|l6f^ 
cases 
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two dimensions |ll|,|l4|], and infinite 
Even in these relatively simple 



rows a rich phase diagram, which is 
not completely elucidated. 

The situation is far from understood in the more rele- 
vant case of three dimensions. There is a general consen- 
sus that in the limit of infinite Hund's coupling (see next 
section) the system has a ferromagnetic phase, although 
phase separation near the Curie temperature at low dop- 
ings has also been found [ p~2|j20|] . If the Hund's coupling- 
is not infinite, or if there are antiferromagnetic superex- 
change interactions present, the situation becomes more 
complicated. There is evidence that an antiferromagnetic 
coupling can induce a first order phase transition at rel- 
atively large dopings pl| . Emphasis on the presence of 
first-order transitions at low temperatures has also been 
done in Refs. p2| , p3[ . In three dimensions, Berry's phases 
in the wave function of the electrons may arise p4] , po| , 
leading to another source of complexity. 

The present work presents a numerical study of the 
double exchange model in three dimensions, using a 
method discussed elsewhere [E6|. We show that the 



interplay between antiferromagnetic interactions and a 
three dimensional structure leads to a rich phase dia- 
gram, where, in addition to ferro- and antiferromagnetic 
phases, other ordered phases with cubic symmetry are 
possible. The model is discussed in the next section. 
Then, we present the calculation at zero temperature 
and the finite temperature results. In the last section, 
we discuss the main features of our findings, and their 
experimental implications. 



II. THE MODEL. 

A. The Hamiltonian. 

We study the double exchange model. In its simplest 
form, electrons in a cubic lattice are coupled to localized 
spins defined at the same lattice sites. The spin of the 
electrons is constrained to be parallel to that of the lo- 
calized spins. This restriction leads to a modulation of 
the electron hopping between lattice sites, which depends 
on the relative orientation of the core spins. In addition, 
one can define a direct coupling between spins at different 
sites. The Hamiltonian is: 



H = — JafJ^ Si 



(<0i&IM. 



■j/ c i c i 



(1) 



where (9i4>i\0j(f)j) is the overlap between the spinors de- 
fined by the polar angle 0, and the azimuthal angle (f>: 

(0i&IMi> = cos | cos | + sin | sin |e-^-*'> . (2) 

We describe the spins S by classical variables, normalized 
to one, \S\ = 1. The system depends on two dimension- 
less parameters only, the filling of the electron band, x, 
and the ratio JafA (notice that our Jaf is opposite in 
sign to the convention of |Q) 

The double exchange Hamiltonian, Eq. (|l]) , can be de- 
rived from the model which describes a lattice of atoms 
with strong intra-atomic Hund's coupling between elec- 
trons in different orbitals. The simplest Hamiltonian 
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which includes this effect is that of the ferromagnetic 
Kondo lattice, 

Ti-FK = ( 4,s C 3,s +h.C.^-J H ^ c\ s Ci^(T ss ,-Si , 

i,j.s i,s,s' 

(3) 

where Jh > 0, and the cr's are Pauli matrices. When 
Jh 2> t, the spin sub-band with an antiparallel ori- 
entation to the core spin lies at high energy, and can 
be projected out. To first order in 1/Jh, one finds an 
antiferromagnetic interaction between nearest neighbor 
core spins, leading to the double exchange Hamiltonian, 
Eq. (Q), with | Jaf| = £ 2 /(4Jh) (note that we have nor- 
malized the core spins). This effect adds to an antifer- 
romagnetic superexchange interaction between the Mn 
t g core spins. Therefore, the most phenomenologically 
interesting range for | J AF |A is 0.008 — 0.15 pi. 



B. Other interactions. 

The model discussed above is supposed to be a reason- 
able starting point to the study of the magnetic prop- 
erties of doped manganites, Lai-^yl^MnOs, where A 
stands for a divalent cation. The electrons of the con- 
duction band arise from the Mn e g orbitals, while the 
core spin is built up from three electrons in the Mn t g 
orbitals. In a cubic lattice, the e g orbitals are doubly de- 
generate. This degeneracy is ignored in Eq. ([!]). Thus, we 
cannot study effects associated to orbital ordering [ p7[ . 
These effects are probably important in explaining the 
different phases with charge ordering observed in these 
materials. The model studied here can be considered a 
first approximation to the study of the magnetic proper- 
ties, in situations where there is no charge ordering and 
the lattice symmetry is cubic. 

The existence of double degeneracy in the e g orbitals 
implies the possibility of lattice Jahn- Teller distortions, 
which have been considered the cause of the breakdown 
of cubic symmetry at low dopings. In addition, dynamic 
Jahn- Teller fluctuations can play a role in some proper- 
ties of doped compounds with cubic symmetry [28 




will ignore these interactions. Previous studies 
gest that the model used here suffices to understand the 
main features of the magnetic transitions in the man- 
ganites, in the doping range 0.1 < x < 0.4 where most 
experiments are done. 

The treatment of the core spins as classical variables 
will overestimate the tendency of the system toward long 
range order. Previous studies Jl(| suggest that this ef- 
fect is small, as results obtained for the physical value 
5 = 3/2 are very close to those found in the limit S — > oo. 
In addition, the only transitions at zero temperature 
found in our work are discontinuous (see next section). 
Hence, quantum spin fluctuations are bounded, and can- 
not change qualitatively the results. 



III. RESULTS. 



A. Phase diagram at zero temperature. 



The phase diagram at zero temperature can be ana- 
lyzed by minimizing the total energy, Eq. ([!]) , as function 
of the spin configuration, for the different values of x and 
the ratio Jaf/£- However, previous studies in one, two 
and infinite dimensions ]s|-|T4|] have shown that the sys- 
tem has a strong tendency toward phase-separation. This 
enforces us to use the Maxwell construction to correctly 
find the phase boundaries. The easiest way of perform- 
ing the Maxwell construction consist in minimizing the 
grand-canonical energy 



n 



GC 



(4) 



as a function of the spin configuration for the different 
values of [i and Jaf- In this way, we obtain the band 
filling, x, as a function of the chemical potential /i/t and 
of Jaf- At the phase boundaries x is a discontinuous 
(growing) function of \x: the jump at the discontinuity is 
the phase-separation compositional range. 

For the above mentioned minimization we will limit 
ourselves to the possible spin textures observed in Monte 
Carlo calculations at finite temperatures (see next sub- 
section) and to those which seem plausible on physical 
grounds. Some of the most stable phases analyzed are 
described in Table I. We have studied only the range 
x < 0.5, due to the particle-hole symmetry of the Hamil- 
tonian (|I|). 

We find the ferromagnetic (FM) and antiferromagnetic 
(G-AFM) phases. The A-AFM configuration is ferro- 
magnetic within planes, and antiferromagnetic between 
neighboring planes. The C-AFM configuration is fer- 
romagnetic along lines, and antiferromagnetic between 
neighboring lines. These phases have already been dis- 
cussed in the literature. The twisted phase interpolates 
between the A-AFM and the C-AFM phases. Island 
phases have been reported in calculations in one and 
two dimensional systems Jl^JlJ]. These are structures 
of spins aligned in one direction with the sense varying 
periodically. For the values of Jaf A considered here, we 
only find the ( ^ , § , 7r) (that corresponds 2x2x1 blocks of 
equal spins) and (^, 7T, it) (3x1x1 blocks), although we 
have tried different combinations 15 1. In these phases, 
the electrons are localized, and the electronic density of 
states is built up of delta functions. Helix phases are or- 
dinary spin density waves. It is interesting to note that 
helix phases are completely irrelevant in the thermody- 
namic limit, but they are remarkably stable on finite lat- 
tices, and should be considered when analyzing numerical 
simulations. 

We find two new phases, intrinsic to the cubic lattice, 
labeled flux and skyrmion in Table I. The spin configu- 
ration is not coplanar, and the spins are parallel to the 
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diagonals of the unit cube (note that the energy is in- 
variant under a global rotation, so that only the relative 
angles are relevant). In the skyrmion phase, the spin 
directions around a unit cube point out from the cen- 
ter, with a hedgehog shape. The induced Berry phase 
can be thought of as generated by an effective magnetic 
monopolc inside each of the unit cubes which build up the 
three dimensional lattice. The charge of the monopoles 
change sign in neighboring cubes, leading to zero total 
charge. The flux within each cube is isotropic, and equal 
to 7r/3 (skyrmion configuration) and 27r/3 (flux config- 
uration) per plaquette. In both cases, the lattice unit 
contains eight sites, and the global symmetry remains 
cubic. 

The electronic dispersion relation for all these phases 
can be calculated analytically, and it is given in Table I. 
The corresponding density of states is shown in Fig. [I] . 
The dispersion relation in the flux and skyrmion phases 
depend linearly on k near k = (±?,±5,±S), and re- 
sembles the dispersion of Dirac massless particles in the 
lattice. 




E 



FIG. 1. Density of states for several spin configurations 

The flux and twisted phases show similarities to the 
pha se with spins at right angles found in two dimensions 

In the twisted phase, for a = b = l/v2, the planes xz 
and yz have the same structure as the 2D pattern studied 
in p2 , pQ , H . The flux phase, on the other hand, has a 
dispersion relation which is the natural extension to 3D 
of the one for the 2D flux phase. Note, however, that 
the flux phase in 2D is coplanar, and all phases in the 
hopping elements can be rendered real by an appropriate 
gauge transformation, while this does not happen in the 
3D phase studied here. 

The stability of the flux and skyrmion phases arise 
from the canting of the spin orientations, coupled to a 
shift in the electronic density of states toward the band 
edges. It is interesting to note that the electronic energy 
of the flux phase is close to that of the C-AFM one, and 



the energy of the skyrmion phase is also close to that 
of the A-AFM phase, for the entire range of electronic 
concentrations . 

The phase diagram at zero temperature is shown in 
Fig. 0. All transitions are first order. The most inter- 
esting result is the stabilization of the skyrmion phase 
for a range of dopings and values of JafA where the 
model is applicable to real materials. The twisted phase 
changes continuously from a 2 = 0.45 (right) to a 2 = 0.77 
(left). Notice that although the A-AFM phase is beside 
the twisted phase, there is a discontinuity in the a pa- 
rameter (a = for A-AFM) . We find island phases (l4| in 
a more restricted range of parameters than in two dimen- 
sions, possibly due to the competition with the flux and 
skyrmion phases. We have not found any tridimensional 
island as (■§,■§,§)■ For more negative values of JafA we 
also find a (§,7T, 7t) island phase starting at —0.19. 



X 




-0.15 -0.10 -0.05 0.00 



FIG. 2. Phase diagram at T — 0. Upper part: chemical 
potential (fi/t) vs. JafA- Lower part: concentration (a;) vs. 
JafA- The results shown in this plot have been obtained in 
a 2400 3 lattice. We have checked that the variations due to 
finite-size effects are smaller than the line widths. 
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Type 



spin direction 



Ferromagnetic (FM) 
A-type (A-AFM) 
C-type (C-AFM) 
G-type (G-AFM) 
Twisted 

Flux 

Skyrmion 
Helix 

Island (^, 7T, 7r) 
Island (^, 7T, 7r) 



Island (§, f ,tt) 



(0,0,1) 
(0,0, (-1) 2 ) 
(0,0, 
(0,0, (-lf+sH-*) 

(a(-l)^,6(-l) z ,0) 
a 2 + b 2 = 1 

((-l)^,(-l)*+*,(-l)^)/V3 

a-ir, (-i)», (_i)»)/v3 

(cos(wz), sin(wz), 0) 

(o,o,(-i)[f]+^ z ) 

(0,0,(-l)[*l+^) 

(o,o,(-i)[*]+^]+ 2 ) 



- 2 E M =iCosfc M 
- 2 E^=iCosfc A , 

— 2 COS &3 



±2 -y/a 2 cos 2 fc3 + 6 2 (cos fci ± cos &2 



±2 y | y X)^=i cos2 fc P ± y 3 cos2 cos2 

—2 X)u=i cos — 2 cos(lj/2) cos(fc 3 +w/2) 

-1,1 
-V2,0,V2 
-2,0,0,2 



Table I. Spin configurations and electronic dispersion relations of the different phases considered in the text. The 
notation [•] stands for the integer part. For a lattice of linear dimension L, ui and k^ can be written as ^-n, with n 
integer (in twisted, flux and skyrmion configurations n is even as the unitary cell is a 2 3 cube). 



B. Finite temperature results. 

We have extended the previous studies to finite tem- 
peratures by using the hybrid Monte Carlo method re- 
ported elsewhere [j26| , which allows to thermalize 16 x 16 x 
16 clusters in reasonable computer time. We have studied 
the model ([!]), on the cubic lattice of side L with periodic 
boundary conditions, using the grand-canonical ensem- 
ble (|). The use of an efficient algorithm has been crucial, 
because around two hundred points of the (T, fi, Jaf/£) 
phase-diagram have been studied. 

Before presenting our results it is worth to comment on 
the strong finite-size effects of this model. Fortunately, 
one can have a good analytical command on them (at 
least at zero t empe rature), by repeating the minimization 
of subsection III A on a finite lattice (this is straightfor- 
ward using Table I). In this way one discovers that the 
phase boundaries change considerably with the lattice 
size. Even worse, for lattice sizes L = 4 and 8, the helix 
phase of minimal frequency to — 2n/L (see Table I) is 
more stable than the ferromagnetic phase in important 
regions of the phase diagram. Also the C-AFM phase is 
stable in the L = 8, although it is unstable in the L — > oo 
limit. On the other hand, the zero temperature phase di- 



agrams of the L = 6 and L = 10 lattices are rather close 
to the infinite volume limit. The only caveat of this lat- 
tice sizes is that they cannot accommodate island-phase 
whose unit cell is not commensurate with them (in par- 
ticular, the (§ , § , 7r) island phase, that has a 4 x 4 x 2 unit 
cell). We have therefore chosen to perform most of the 
numerical simulations on the 6 3 cluster and, only when 
necessary for elucidation, we have employed as a larger 
cluster the 10 3 . Exception to this rule has been the sim- 
ulation of the 8 3 in the range -0.29 < JafA < -0.18 
range, to check for the presence of the 2x2x1 island 
phase. 

In the simulation, we have used the perfect action for- 
mulation |2q| , with A = 0.25. The molecular dynamics 
parameters have been chosen to ensure an acceptance 
rate greater than an 80%. The selection depends mainly 
on the temporal lattice size (inverse temperature). Typ- 
ical values used are N T = 20, At = 0.02 for T/t = 1/20, 
and N T = 50, At = 0.01 for T/t = 1/45. We have sys- 
tematically controlled the correct thermalization of our 
data, by comparing the results of hot and cold starts. 
In some cases, in particular at T = t/100, we have also 
compared the results of ordered starts with a slow an- 
nealing from high temperatures. The total run-length 
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(up to 2000 trajectories in some points) has always been 
at least larger than three times the thermalization time. 
The typical size of our statistical errors can be found in 
figure 0. The total CPU time devoted to this project has 
been the equivalent of 10 years of a Pentium III at 800 
MHz. 



in Table I, it follows that the appropriate pseudo order- 
parameters are: 



0.8 



0.6 



0.4 



0.2 



♦ o a 

A G] 
a X 



-o- .-o- -o- -o — o- 



-0.15 



-0.05 



-0^— J — O 1 — O — -O 



0.00 



-0.10 

J„/t 

FIG. 3. Monte Carlo measures of several quantities at 
T = t/45 for fi = (upper side) and = —1.5, on the L = 6 
lattice. 

As we have said in the introduction, our scope has been 
the determination of the magnetic phase diagram of the 
model. Since one cannot break a symmetry on a finite 
lattice, it is necessary to use pseudo order-parameters for 
obtaining the phase diagram p2fl . We have defined our 
pseudo order-parameters in terms of the spin structure 
factor (V = L 3 is the total number of spins), 



(5) 
(6) 



In the above equation, (•) stands for the thermal average. 
In a paramagnetic phase, G(fc) will be of order 1/V for 
all fc. On the other hand, on an ordered phase it will 
be of order one for a very small set of fc's (see below), 
and of order 1/V for the rest of wave numbers. If one 
now calculates the structure factor for the configurations 



G(fc) = (|S(fc)' 2 " 



G r 


= 0(0,0,0), 






(7) 


G s 


= G(7T, 7T,7t), 






(8) 


G\ 


= G(7T,7T,0) ■+ 


-G(7T,0,7T) + 


G(0,7r,7r), 


(9) 


G p 


= G(7T,0,0)-t 


G(0,7T,0) + 


G(0,0,tt), 


(10) 


G'h 


= G(f ,0,0) 


+ G(0,^,0) 


+ G(0,0,f ). 


(11) 



Let us explain them in details. In a ferromagnetic phase, 
only Gf is of order 1, while for a G-AFM phase, the 
appropriate pseudo order-parameter would be G s . For 
an A-AFM phase, the selected order parameter will be 
G p (we sum for fc = (tt, 0, 0), (0, tt, 0) and (0,0, tt) be- 
cause we do not know a priori the direction of the AFM 
lattice planes). A skyrmion phase will also have G p 
as pseudo order-parameter, and so the question arise of 
how to differentiate a skyrmion from an A-AFM phase. 
The easiest way to do that is to control that for the 
skyrmion phase the three quantities G(ir, 0, 0), G(0, tt, 0) 
and G(0, 0, 7r) are simultaneously large (we have also nu- 
merically checked the orthogonality of the Fourier trans- 
formed spin field at the three fc values). Also the C-AFM 
phase and the flux phase have a common pseudo order- 
parameter, namely G\, but one can differentiate both 
phases in exactly the same way. A twisted phase will 
have as pseudo order-parameter both G p and G\. One can 
distinguish the twisted phase from a tunneling phenom- 
ena between an A-AFM and a C-AFM phases, checking 
again that both pseudo order-parameters are simultane- 
ously large (we have also checked the orthogonality of 
the Fourier transformed spin-field at the two fc for which 
the structure factor is not negligible). Finally the helix 
phase is signaled by the Gh order parameter, while the 
structure factor for the (k±, k2,ks) island phase is heavily 
peaked at fc = (fci, fe, £3). 

The results for the G a pseudo order-parameters and 
for the equilibrium fermionic density at T/t = 1/45 are 
shown in Fig. for [i = (half filling, x = 0.5) and fi = 
— 1.5t which includes the interesting region of x ~ 1/4 — 
1/3. The numerical results shown in this plot allow us to 
identify the FM, G-AFM, twisted and flux phases, and 
the onset of the skyrmionic phase. 

For /.i = 0, the flux-twisted transition is clearly dis- 
continuous. In addition, there is a clear metastable be- 
havior (the duplicated points correspond respectively to 
two metastable phases) . At the right end of the twisted 
phase, G p and G\ vanish abruptly and there are clear 
signs of metastability with a phase where Gf is small but 
non- vanishing. To check if it is a finite size effect, we 
have performed some simulations in this region with a 
larger lattice (10 3 ), obtaining that Gf is much smaller 
(and G p and G\ vanish). Thus we conclude that there is 
a first order transition at this temperature between the 
twisted and the paramagnetic phases. On the contrary, 
the ferromagnetic-paramagnetic transition is smoother. 
We recall that in Ref. p6[ it is shown that at Jaf = 0, 
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T/t = 0.14 the transition is second order. Although we 
find a strengthening of this transition when | Jaf| grows, 
to confirm that it actually becomes of the first-order 
would require a careful finite-size scaling analysis that 
is beyond the scope of this work. 



For /i = —1.5 (see Fig. ||) we find a very different 
landscape. At the left part we find an almost saturated 
G-AFM region (notice that x as and the fermions es- 
sentially play no role). In the center, we find a region 
with x ~ 0.25 where G p is about 0.05 and the others are 
much smaller. We interpret this result as the onset of the 
skyrmionic region (that appears clearly at lower temper- 
atures, see Fig. ||). The start of the ferromagnetic region 
is much abrupter than at fj, = 0, but as said above, to 
conclude that it is first order, a study in different lattice 
sizes is mandatory. 
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FIG. 4. Phase diagram obtained with Monte Carlo mea- 
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of the circles are proportional to the associated G a - For clar- 
ity at each point we plot only the maximum G a except in the 
cases where G p /G\ £ [0.25, 1], that correspond to the twisted 
phase. 
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Results for the different G a order parameter for fx = 
— 1.5f are shown in Fig. ^[ Skyrmion correlations are 
suppressed at relatively low temperatures, but for T = 
i/100 we find a clear skyrmionic structure (G p = 0.7 
— 0.8, with the appropriate orthogonality of the Fourier 
components). 
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FIG. 6. Phase diagram at T = t/45. The solid lines corre- 
spond to T = for a L = 6 lattice. 



To explore the phase diagram, we have simulated at 
several values of T. The phase diagram at half filling is 
plotted in Fig. ||. The paramagnetic region separates the 
ferromagnetic from the twisted phase for T/t > 1/45 but 
it is difficult to know the point where it finishes (maybe 
T = 0) . As the ferromagnetic magnetization goes to zero 
rather smoothly, to obtain a more precise value of the 
transition point we have used data from larger lattices 
(mainly L — 10). The dashed lines have been draw ap- 
proximately at the phase transitions. 



We have combined all our results at T = t/45 in the 
phase diagram shown in Fig. |6[ The results are con- 
sistent with the zero temperature calculation for the 6 3 
lattice, whose transitions are plotted as continuous lines. 
The FM, G-AFM and Flux phases (and the twisted one 
near half filing) are less affected by finite temperature 
effects than the other phases, that become paramagnetic 
for temperatures as small as T = i/45. Notice also that 
in the 6 3 lattice the island phases are so narrow that they 
are hardly observed in a simulation. However, we have 
checked that in a 8 3 lattice at \i = 0, T — t/45 the sys- 
tem is in the 2x1x1 phase for JafA <= (—0.29, —0.18). 
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We have also simulated a 8 3 lattice in several values of 
(fi, Jaf) that at T — correspond to the 2x2x1 region, 
observing that even at T = t/45 the system is disordered. 



C. The skyrmion phase. 

As mentioned above, we find a novel phase, the 
skyrmion phase, for a realistic range of dopings and anti- 
ferromagnetic couplings. We now discuss some physical 
features of this phase. Similar properties can be expected 
in the flux phase. 

i) The symmetry of the magnetic phase is cubic. This 
phase is compatible with experiments done in the range 
x ~ 1/4 — 1/3 which show phases with cubic structure 
but no macroscopic magnetization |[33| . 

ii) The magnetic arrangement is, at first sight, similar 
to that expected in a canted phase. Measurements of the 
microscopic magnetic structure will probably be indistin- 
guishable from those expected for a canted phase. 

iii) The magnon spectrum will resemble that of other 
systems with non collinear spins [ p^j35| . The dispersion 
relation will be linear at low energies. 

iv) The phase is metallic. However, the chemical po- 
tential, at the appropriate dopings, coincides roughly 
with the position when the density of states has a pro- 
nounced minimum (see Fig. ^ ). Hence, it cannot be ruled 
out that residual disorder will induce insulating behav- 
ior. If this is the case, it will be behave in a similar way 
to the canted insulating phases reported in the literature 
(see, for instance jS6|). 

v) The unusual spin configuration leads to an effective 
coupling between an external magnetic field and a charge 
density wave. An applied field tilts the spins. The sign 
of the change in the angle between neighboring spins de- 
pends on the relative orientation of the spins and the 
magnetic field. A reduction in the angle leads to an in- 
crease of the hopping, and the opposite happens if the 
angle is increased. Hence, electron charge will be dis- 
placed from the weak bonds to the strong ones, and a 
charge density wave will arise. This feature can clearly 
identify the skyrmion phase. 



IV. CONCLUSIONS. 



We have expressed our results in terms of the elec- 
tron hopping, t = W/12, where W ~ 2eV is the width 
of the conduction band, and the effective antiferromag- 
netic coupling, Jaf> arising from superexchange interac- 
tions between the Mn t g orbitals (|Jaf|A ~ 0.005 - 
0.012 |^]) and from the finite value of the intra-atomic 
Hund coupling, Jh . It follows that the most phenomcno- 
logically interesting range for | JafIA i s 0.008 — 0.15 
Using the previous values, we find that the highest tran- 
sition temperature (x = 0.5, Jaf = 0) is close to 350 K, 
and therefore most of the phase transitions reported in 
this work occur in the — 300K range, in agreement with 
experiments. 

It is somewhat surprising that conventional spin den- 
sity wave structures are not stable. The three dimen- 
sional phase diagram shares some features with solutions 
of the double exchange model in one, two and infinite 
dimensions, like phase separation. A number of phases, 
however, have no counterparts in other dimensions. 

We have identified, among others, a phase where the 
spins are locally arranged in a hedgehog manner, labelled 
skyrmion phase. The existence of this structure is prob- 
ably intrinsic to the double exchange mechanism, and 
cannot be realized in spin systems with short range cou- 
plings. The skyrmion phase is stable for dopings x ~ 0.2 
— 0.3 and antiferromagnetic couplings | JafIA ~ 0-05 — 
0.07. These values are realistic for doped manganites in 
the range where Colossal Magnetoresistance effects are 
found. While it is probably difficult to characterize this 
phase experimentally, it will be interesting to verify its 
existence. 
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We have analyzed the phase diagram of the double 
exchange model in a cubic lattice, both at zero and at fi- 
nite temperature. Our results at finite temperature have 
been obtained using an efficient Hybrid Monte Carlo al- 
gorithm |26| l , which has been crucial for the exploration of 
a three dimensional phase diagram. We have found that 
the competition between the ferromagnetic interaction 
mediated by the conduction electrons and antiferromag- 
netic couplings which arise from superexchange effects or 
from the underlying intra-atomic Hund's coupling lead 
to the existence of a variety of phases. 
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